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We study for what specific values of the theoretical parameters the axion can form 



O 



the totality of cold dark matter. We examine the allowed axion parameter region 
in the light of recent data collected by the WMAP5 mission plus baryon acoustic 



O 

£S) ■ oscillations and supernovae, and assume an inflationary scenario and standard cos- 

mology. We also upgrade the treatment of anharmonicities in the axion potential, 
\ which we find important in certain cases. If the Peccei-Quinn symmetry is restored 

Oh; ' J ' 

an axion mass m a = (85 ± 3) /xeV makes the axion 100% of the cold dark matter. 

If the Peccei-Quinn symmetry is broken during inflation, the axion can instead be 

100% of the cold dark matter for m a < 15 meV provided a specific value of the 

(N 

J> ■ initial misalignment angle 0{ is chosen in correspondence to a given value of its mass 

l>. 

t — . m a . Large values of the Peccei-Quinn symmetry breaking scale correspond to small, 

m 

perhaps uncomfortably small, values of the initial misalignment angle 9i 

O , 

Q\ PACS numbers: 14.80.Mz, 95.35.+d 

O ' 

> 

I. INTRODUCTION 

The recent measurements by the WMAP mission [1] have established the relative abun- 
dance of dark and baryonic matter in our Universe with great precision. About 84% of the 
content in the Universe is in the form of cold dark matter (CDM), whose composition is 
yet unknown. One of the most promising hypothetical particles proposed for solving the 
enigma of the dark matter nature is the axion [2] [3]. This particle was first considered in 
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1977 by R. Peccei and H. Quinn [4] in their proposal to solve the strong-CP problem of the 
QCD theory. Although the original PQ axion is by now excluded, other axion models are 
still viable [5-8]. 

The hypothesis that the axion can be the dark matter particle has been studied in various 
papers (see e.g. [9-16] and the reviews in [17, 18]). Here we examine the possibility that 
the invisible axion may account for the totality of the observed CDM, in the light of the 
WMAP5 mission, baryon acoustic oscillations (BAO) and supernovae (SN) data. We also 
upgrade the treatment of anharmonicities in the axion potential, which we find important 
in certain cases. We consider invisible axion models, in which the breaking scale of the PQ 
symmetry f a is well above the electroweak scale. The axion parameter space is described by 
three parameters, the PQ energy scale /„, the Hubble parameter at the end of inflation Hi 
and the axion initial misalignment angle 0j. 



The solution to the strong CP problem proposed by R. Peccei and H. Quinn [4] introduces 
a new U(1)pq symmetry in the theory of strong interactions. The Peccei-Quinn theory in 
its simplest version depends on a unique parameter, the energy scale at which the U(1)pq 
symmetry is broken. The axion field a(x) originates from the U(1)pq symmetry breaking 
at a temperature T ~ f a . Axions are the quanta of the axion field [2, 13]. 

The temperature-dependent axion mass arises through instanton effects and is given by 



We set A = 200 MeV [20], and take the model-dependent factor b = 0.018, consistently with 
previous work [15-17]. The zero-temperature mass m a = m a (T = 0) is [2] 



where z ~ 0.56 and m n and f n are the pion mass and decay constant respectively. The integer 
N represents the U(1)pq color anomaly index; in this paper we set N — 1, consistently with 



II. AXION PROPERTIES 



[19] 




(1) 





[18]. 
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In the following, we set Q 2 = j^f v m n = (78.7 MeV) 2 . From m a {T) in Eq. (1) we see 
that the axion is essentially massless down to a temperature T = O(GeV). This fact is 
extremely important when axions are studied in a cosmological context as in this paper. 
Defining the misalignment angle as 

0{x) = (3) 

J a 

the evolution of the zero mode of the dynamical field 9(x) in a flat Friedmann- Robertson- 
Walker metric is 

+ 3H(T) 9+^^ = 0. (4) 

A dot indicates a derivative with respect to time, H(T) is the Hubble parameter and the 
axion potential is 

V{9)=ml{T)fl{\- cos9). (5) 

For small 9, the potential is approximately harmonic, V{9) ~ \m 2 a {T) f 2 9 2 . In this case the 
equation of motion Eq. (4) becomes 

9 + 3H(T)9 + m 2 a (T)9 = 0. (6) 

When T>A, the axion is massless, see Eq. (1), and a solution to Eq. (6) is 9 = 0j = const, 
where 9i is the initial value of the misalignment field, the so-called initial misalignment angle. 
The axion field is frozen at the value of 9i from the onset of production until the axion mass 
becomes relevant and the axion field begins to oscillate. This occurs below a temperature 
7\ = O(GeV) defined by 

3^(T 1 )=m (T 1 ). (7) 
To find Ti in terms of m a , we use the Friedmann equation for a radiation-dominated Universe, 

"8tt 3 ^(T)" 



H(T) 



90MJ, Z 



1/2 rp2 



T 2 » l.m^JT) (8) 

Mpi 



where g*(T) is the total number of effective degrees of freedom at temperature T [20]. For 
the range of temperatures of interest in this paper, one has 



9*(T) 



61.75, for T > A, 

10.75, for A > T > 4MeV, (9) 
3.36, for T < 4MeV. 
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Inserting Eqs. (1) and (8) in Eq. (7) gives 

( im -^ |A4 ) 1/6 = 618Mevf !0^V /6 , T > A, 

^4.98 v /^m)^ { fa J > ~ ' 



68.1MeV 



( 



10 18 GcV 

/a 



(10) 



T< A. 



^ V 4.98-^/9,^1)' 

Since we expect / a ~ 10 12 GeV, the axion starts oscillating at T x ~ 618 MeV. However, we 
will see that it is also possible for f a to be of order the GUT scale (~ 10 16 GeV) or even the 
Planck scale (~ 10 19 GeV) and coherent axion oscillations start later. 



III. WMAP BOUNDS ON ifj 

The Hubble expansion rate at the end of inflation Hi can be constrained using data from 
WMAP5 plus BAO and SN. The curvature perturbation spectrum A^(A; ) at fixed wave 
number k = 0.002Mpc _1 has been measured as [1] 

A^(k ) = (2.445 ± 0.096) x 10~ 9 . (11) 

The tensor-to-scalar ratio r has been constrained to be 

r = ^f^4 < 0.22 at95%CL. (12) 

Eqs. (11) and (12) can be combined to give an upper bound on the spectrum of primordial 
gravitational waves of about 

A 2 h (k ) < 5.38 x 10~ 10 . (13) 

Expressing A^(A; ) in terms of Hi, 

A^„) = ™jr, (14) 

leads to an upper bound on Hi, 

Hi < 6.29 x 10 14 GeV. (15) 

A lower limit on Hi comes from requiring the Universe to be radiation-dominated at 
T ~ 4 MeV, so that primordial nucleosynthesis can take place [21]. Equating the high- 
est temperature of the radiation 

Tmax ~ {T 2 RH HiM Pl fl\ (16) 
to the smallest allowed reheating temperature T RH = 4 MeV gives 

Hi > H{T RH ) = 7.2 x 10~ 24 GeV. (17) 



5 



IV. BOUNDS ON AXION FLUCTUATIONS 

If the axion energy scale is lower than the Gibbons-Hawking temperature T GH = Hj/2n 
[22], 

fa < T GH , (18) 

the PQ symmetry is restored after inflation (Scenario I). When the Universe expands and 
cools down to a temperature ~ f a , the PQ symmetry breaks again. Different values of 
the misalignment angle 9i(x) are present within one Hubble volume [20], giving rise to 
fluctuations that are adiabatic as observed in the CMB spectrum. 

On the contrary, non-adiabatic fluctuations are generated when (a) the PQ symmetry 
breaks before the end of inflation and (b) it is not restored afterwards (Scenario II). Condition 
(a) requires 

fa > T GH . (19) 

Condition (b) requires 

fa > T MAX , (20) 
where T M ax is given in Eq. (16). Thus, non-adiabatic axion fluctuations arise when 

f a > m&x{T GH ,T MAX }. (21) 

Theories where the inflaton decays into fermions favor T RH ^ 10 12 GeV [23], which is 
consistent with assuming T M ax < T GH . In this scenario Eq. (19) suffices. The initial 
misalignment angle 6i has a single value within one Hubble volume, since it was causally 
connected at the onset of production. The axion mass is negligible, see Eq. (1). Quantum 
fluctuations Sa(x) in the axion scalar field a(x) have variance [22] 

<IM*)| 2 > = (§) 2 - (22) 



It follows that fluctuations in the misalignment angle field 9(x) = a(x)/ f a have variance 

<4 =(&)'■ (23) 



>/«, 

Under condition Eq. (21) axion isocurvature perturbations are present during inflation 
and are constrained by WMAP5. Defining the power spectrum of axion perturbations 
A aO) = (\$Pa/pa\ 2 ), one finds 



= ^hr (24) 

u i J a 
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The axion entropy-to-curvature perturbation ratio is then 

or, introducing the axion adiabaticity ato(ko), 

Ag(fcp) = a (fco) 
A|(fc ) l-ao(fco)' 
The adiabaticity a is constrained by WMAP 5-year data to be 



(25) 



(26) 



a < 0.072 at 95% CL. (27) 
Using the value for A^(/c ) in Eq. (11) this bound can be rephrased as 

< 4.17 x 10~ 5 . (28) 

ViJa 

V. PRESENT AXION ENERGY DENSITY 

The major interest for axions in astrophysics is that it is possible for these particles to 
account alone for all of the observed CDM. If this is true, axions must be in highly non- 
thermal equilibrium and probably form a Bose-Einstein condensate [18, 20]. The leading 
mechanism for producing such an axion population is the misalignment production (see 
e.g. [16-18, 20] and references therein). Another contribution that is important in Scenario 
I is axion production from string decay [18]. 

According to the misalignment mechanism, the axion number density at temperature Ti 
is given by [18, 20] 

C(H) = \m a {T,)mm)) X . (29) 
Here the factor x models the temperature dependence of the axion mass around Ti and 
depends on the number of quark flavors Nf that are relativistic at Ti [24]. We take x — 
1.44, consistent with Nf = 3. The function f(0i) accounts for anharmonicity in the axion 
potential, i.e. for a solution to the full axion field Eq. (4) instead of Eq. (6). The function 
f(9i) is of order one for \6i\ ^3 and is logarithmical divergent for 6i — > n. It is discussed at 
the end of this Section. 

The number density at present time is found by imposing the conservation of the comoving 
axion number density in the form [20]: 

'n™ s (T)' 



s(T) 



0, (30) 
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where s(T) the entropy density. One has 

s(T) = %9*s{T)T\ (31) 

We use the approximation of the same number of degrees of freedom for the entropy and 
the total energy density for all g*s(T) with T > 4MeV; the present value for <?*s(T ) = 3.91 
differs from g*(T ) = 3.36 [20]. Thanks to Eqs. (29), (30) and (31), the present axion energy 
density p™ s (T ) = m a n™ is {T ) is 

Pa lS (To) = ma ^ ( g i ) ) 8(r ° ) / a 2 (g?/(gO>X- (32) 

Dividing the last equation by the critical density p c = ^HlMpJ^n and using Eq. (1) for 
m a (Ti), the cosmologically relevant ratio f2™ s = p a j ' p c is 

nr = P^<$m))x 1 6( * )4 ' /a f L (33) 

2 ^ Tl ) [ 1, fa>fa- 

Or, if we insert the expression for Ti previously computed in Eq. (10) and substitute the 
numerical values x — 1-44, b = 0.018, we obtain 

arh2 = WmqmMmk&r'*, fa z h, (34) 

|o.0051(^ 2 / (^))( T ^)3/2 ? fa >f a 

The scale /„ is computed by equating the two expressions for Vt" as h 2 above: 

f a = 9.91 x 10 16 GeV. (35) 

The two cases in Eq. (34) reflect the dependence of the axion mass m a (T) on the temperature 
T in Eq. (1). 

The factor (9ff(9i)) assumes different values in the two cases in which the axion is formed 
after inflation, f a < T GH (Scenario I), or is present during inflation, f a > T GH (Scenario II). 
In Scenario I, the variance is zero because there are no axion quantum fluctuations from 
inflation, but 9i is not uniform over one Hubble volume, so 6f is averaged over its possible 
values as 

@m)) = ^ fjumde,. (36) 

If f(0i) is assumed to be 1, one finds the usual result (6f) = 7r 2 /3. If f(9i) is taken to be 
the more realistic expression in Eq. (40) below, one finds 

@m)) = 8.77. (37) 
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Other expressions for f(9i) lead to similar values [24]. 

In Scenario II, the initial misalignment field has a single value 9i over one Hubble volume, 
so = 9i. The variance of the misalignment field is given by Eq. (23), so 

2" 



el 



Hi 
2nf a 



f(0i). (38) 



The initial misalignment angle 9i in this scenario can assume any value between — n and ir. 
When 6i <C 1, the axion energy scale f a ^> 10 12 GeV if the axion is to be the CDM particle 
[25]. 

We now discuss the anharmonicity factor f(9i). This factor accounts for the corrections 
that must be made to the solutions of the linear Eq. (6) to obtain the solutions of the full 
non-linear Eq. (4). The anharmonicity factor f(9i) has the following limiting behavior: 

hm m) = l; lim f(9 l ) = +oo. (39) 

The problem of finding a specific shape for /(#«) has been considered by many authors. 
Turner [24] integrates Eq. (4) numerically and describes how the anharmonicity factor can 
be computed, but does not give an explicit formula. Lyth [26] follows the idea in [24] and 
performs an explicit calculation, obtaining the behavior f(9) ~ ln L175 (l — ^/7r) for 9i > 0.97r; 
he also comments that his result differs from Turner's by a factor of two. The exponent 1.175 
comes from Lyth's dependence of m a (T) on T, which is m a (T) ~ T" 3 7 , different from ours. 
Finally, Strobl and Weiler [27] and Bae et al. [28] perform a more precise numerical analysis 
following [24], and confirm the result in [26] for the behavior of f(9i) around 9i = it. 

In this paper, we adopt the following analytic expression for /(#«), which extends Lyth's 
formula to values of 9i < 0.97T and symmetrically to negative values of Of 



In 



7/6 



(40) 



This expression has the limiting behaviors in Eq. (39) and is analytic in the range (— 7r,7r), 
so the computation of the level curves in Fig. 1 below can be carried out analytically. The 
different exponent in the power of the logarithm, 7/6 ~ 1.167 instead of Lyth's 1.175, comes 
from the different dependence of the axion mass m a (T) on T in Eq. (1), m a (T) ~ x~ 3 ' 7 in 
Lyth's paper and m a {T) ~ T -4 in this paper. 

Other authors find different expressions for the axion energy density. 
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Hertzberg et al. [16] use different values for x, which range from x — 1 (which they call 
a "moderate" value) to x = 1/20 (which they call a "conservative" value). However, the 
authors in [16] do not account for anharmonicities in the axion potential near 0i = tt. As 
shown in Section VI, the anharmonicity factor is essential to obtain the correct behavior 
of the isocurvature fluctuation bound at relatively small values of f a , because the initial 
misalignment angles that give the correct axion density are not very small. We choose 
X = 1.44, but most importantly we differ from [16] in that we account for the behavior of 
the anharmonicity function /(#»). With this prescription, it turns out that the figures shown 
in [16] are modified when the function ceases to be of order one. This happens for 

f a ^ 10 11 GeV, as shown in Fig.l below. 

Sikivie [18] studies the axion energy density in the case T\ > A, and finds Q a = 
0.15 6f (/ a /10 12 GeV) 7 / 6 . His expression differs from ours in Eq. (34) because of the different 
numerical factors used in [18], namely x =l,b = 5/12 and T\ = lGeV(10 12 GeV// (l ) 1 / 6 . 

VI. RESULTS 

In this section we show the region of axion parameter space in which the axion is 100% 
of the cold dark matter. In other words, we assume 

n a h 2 = n CDM h 2 = 0.1131 ±0.0034 at68%CL. (41) 

Here VL C dm is the density of cold dark matter in units of the critical density and h is the 
Hubble constant in units of 100 km s _1 Mpc _1 . 

The parameter space is labeled by the axion energy scale f a , the initial misalignment 
angle 9, h and the Hubble parameter during inflation Hj. Results are shown in Fig. 1. 

The region labeled "Tensor modes" shows the constraint on Hi in Eq. (15) coming from 
the WMAP5 plus BAO and SN observations [1]. The newly launched PLANCK satellite 
will improve the actual measurement on the tensor-to-scalar ratio r by at least one order of 
magnitude [29]. If PLANCK will not detect gravitational waves and r < 0.02, Hi will be 
bound by Hi < 2 x 10 14 GeV. This forecast measurement is shown by a vertical dashed line 
labeled "PLANCK". 

The region labeled "White dwarfs cooling time" is excluded because there one would have 
an excessively small cooling time in white dwarfs. These methods set a limit on the axion 
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H r [GeV] 

FIG. 1: Region of axion parameter space where the axion is 100% of the cold dark matter. The 
axion mass scale on the right corresponds to Eq. (2) with U(1)pq color anomaly N = 1. When 
the PQ symmetry breaks after inflation (f a < Hi/2ir), the axion is the CDM particle if f a = 
(7.27 ± 0.25) x 10 10 GeV, or m a = (85 ± 3) fJ,eV, which is the narrow horizontal window shown 
on the right (we plot a 3a window to make it visible). If the axion is present during inflation 
(f a > Hj/2ir), axion isocurvature perturbations constrain the parameter space to the region on 
the top left, which is marked by the values of 9i necessary to obtain 100% of the CDM density. 
Other bounds indicated in the figure come from astrophysical observations of white dwarfs cooling 
times and the non-observation of tensor modes in the Cosmic Microwave Background fluctuations. 
Dashed lines and arrows indicate the future reach of the PLANCK satellite and the AD MX and 
CARRACK microwave cavity searches. 

energy scale of [30] 

/ fl >4x 10 8 GeV. (42) 

Assuming N = 1 in Eq. (2), this corresponds to m a < 15 meV. 
The line 

U = T GH = Hi/2n (43) 
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divides the region where the PQ symmetry breaks after inflation, f a < T GH (Scenario 1), 
from the region where the axion field is present during inflation, f a > T GH (Scenario II). 

The region marked as ADMXI has been excluded by a direct search of axions CDM in 
a Sikivie microwave cavity detector [31] by the ADMX Phase I experiment [32, 33]. The 
window shown corresponds to 1.9/xeV < m a < 3.3 /xeV, valid for the KSVZ axion model. 
A narrower DFSV axionic window, not shown in the figure, has also been ruled out. The 
dashed line labeled "ADMXII" shows the forecast axionic region to be probed in the ADMX 
Phase II, which would search for axions with mass up to 10/ieV, or f a = 6.2 x 10 11 GeV. The 
proposed CARRACK II experiment is a cavity search that will look for axions with mass 
up to 50 fieV [34]. 

For Scenario I, since (0ff(6i)) does not depend on Hi or the value of f2™ s depends on 
f a only. From Eq. (34) with f a < f a and Eq. (37), we find 

nrh2 = 2 - 07 {iJhv) 7/6 - (44) 

As pointed out in [18] , axionic string decays are an important mechanism of axion production 
in Scenario I (in Scenario II all string defects are washed out by inflation). The present 
density of axions produced in string decays is [18] 

n "' fc2 = °- 34 (lo4v) 7/6 - (45) 

The total axion energy density is therefore 

n a h 2 = nr s h 2 + nfh 2 = 2.41 (y^ri^v) ^ • ( 46 ) 

The value of f a such that the axion is 100% of the cold dark matter, fl a h 2 = flcDMh 2 , is 
then 

f a = (7.27 ± 0.25) x 10 10 GeV. (47) 

This band is drawn in Fig. 1 as the horizontal window on the lower right. The result in 
Eq. (47) is one order of magnitude lower than the usually quoted value for the PQ energy 
scale, fa ~ 10 12 GeV. This difference by one order of magnitude comes from equating Q a with 
the precise result in Eq. (41) for Qcdm, that is Q Q ~ 0.1, whereas seminal works in axion 
cosmology computed f a from the equation Q a — 1. We also include in Q a the contribution 
from axionic strings decay, Eq. (45), and the careful derivation of Eq. (44). 
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Assuming a U(1)pq color anomaly N — 1, the Q a = VL C dm band corresponds to an axion 
mass 

m a = (85 ± 3) /ieV. (48) 

In Scenario II, axion isocurvature fluctuations are present and lead to the bound in 
Eq. (28). Together with the expression for Q™ is in Eq. (34) and the condition of 100% 
CDM in Eq. (41), the adiabaticity bound excludes the shaded region in the center of Fig 1. 
The PLANCK satellite is expected to improve the current bounds on the axion isocurvature 
fluctuations by at least one order of magnitude [29]. The dashed line on the left of Fig. 1 
shows the new bound on the allowed region when Eq. (27) is replaced by ao < 7 x 10 -3 . 

The leftmost boundary of this region contains two kinks and can be approximated with 



fa 



7.64 x 10 9 GeV ( w ,% cY ) , < 9.96 x 10 6 GeV, 

1.97 x 10 13 GeV (jp^y) 12/5 , 9.96 x 10 6 GeV < H 1 < 3.43 x 10 9 GeV, (49) 
^6.74 x 10 18 GeV ( 1Ql0 g/ GeV ) 4 , > 3.43 x 10 9 GeV. 



The upper kink occurs at f a = f a and is due to the change in the dependence of the axion 
mass on the temperature, Eq. (1). The lower and smoother kink around f a ~ 10 11 GeV 
arises from the fact that the anharmonicity function f(6i) differs from one at values of f a 
smaller than 10 11 GeV, see Eq. (40). Notice that the simple proportionality f a oc Hi at 
small Hj < 10 7 GeV is independent of the detailed form assumed for the function f{6i) near 
9i = n, and derives in a straightforward way from Eq. (28). 

In the remaining region on the left of Fig. 1, the axion can be 100% of the cold dark 
matter, provided the value of the initial misalignment angle 8{ is chosen appropriately. 
The 9i contours in the figure indicate the appropriate values of 9i for a given Hj and f a . 
Notice that the 9i contours are horizontal, i.e. independent of Hi, since in that region the 
contribution from in Eq. (38) is negligible compared to Of. This allows us to show the 
full relation between f a and ^ imposed by the constraint Q a = VL C dm- We find 



h 2 \ 6 / 7 f < ? a > n nm 

Ja I U.236 2 Wi)) ' Ja ~ Ja UL Ui ~ U ' UU1 > 



10 12 GeV 1 / n h 2 .2/3 > , < 

U.0051 6j f(6i) ) ' Ja ~ Ja Ui V l ^ U.UUl. 



(50) 



This is illustrated in Fig. 2 for Vt CDM h 2 in Eq. (41). 



13 



OS 



1.0 



0.8 



0.6 



0.4 



0.2 



0.0 



m a [eV] 

10" 3 10" 4 10" 5 10" 6 10" 7 10" 8 10" 9 



da < ^CDM 



&a > ^CDM 



10 9 10 10 10 11 10 12 10 13 10 14 10 15 10 16 
fa [GeV] 



FIG. 2: The misalignment angle 9i necessary for the axion to be 100% of the cold dark matter in 
Scenario II (f a > Hi/2tt), as a function of f a - Above the curve, O a > CIcdm- For f a > 10 17 GeV, 
one has 9i ~ 0.001(/ a /10 17 GeV) _3//4 ; in particular, for f a > 10 19 GeV, the initial misalignment 
angle 9% has to assume values 9i ^$ 10~ 5 . 

For f a > 10 17 GeV (m a < lO" 10 eV or B % < 0.001), one has /(0 4 ) ~ 1 and Eq. (50) 
simplifies to 

f \ "3/4 
J a 



Bi ~ 0.84 x 10 



-3 



10 17 GeV 



or 



Bi ~ 1.2 x 10~ 3 (- 



m a \ 3 / 4 



for f a > 10 17 GeV, 



for m a < 10~ 10 eV. 



(51) 



(52) 



10- 10 eV7 

In particular, for f a > 10 19 GeV, the initial misalignment angle Bi has to assume values 
Bi ^$ 10 -5 . This was also noted in [17, 28]. These small values of Bi may be uncomfortable 
in a cosmological scenario. 

In the other limit of Bi ~ tt, the form of the function f(Bi) assumed in Eq. (40) gives 

en 



ir-Bi- —e- c/s % for f a < 2 x 10 10 GeV, 



(53) 



with C = 7.48 x 10 10 GeV. So, as Bi approaches it from below, the corresponding f a ap- 
proaches 0. This gives rise to the linear dependence of /„ on Hj in the lower left corner of 
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Fig. 1. 

VII. DISCUSSION 

We have seen that depending on the ratio f a /Hj, two scenarios are possible for the axions 
to be 100% of the cold dark matter, which we called Scenario I (f a < Hj/2tt) and Scenario II 
(f a > Hi/2n). In Scenario I, the Hubble scale Hi is bounded from below to Hi > 10 11 GeV 
(leftmost edge of the Q a = ^cdm horizontal window in Fig. 1). This window can disappear 
completely if the limit from tensor modes moves to the left beyond Hi > 10 n GeV. If the limit 
on r would become more stringent than r > 10~ 8 , Scenario I would have to be abandoned 
in favor of Scenario II. 

We remark that for very large values of f a in Scenario II, the initial misalignment angle 
9i has to be chosen very small for the axion to be 100% of the CDM, see Eq. (50) and Fig. 2. 
This may undermine the axion field as a dynamical solution to the strong CP problem, in 
that it would have to be fixed to a small value as an initial condition. 

Grand Unification Theory (GUT) models that contain axions predict that f a should be 
of the order of the GUT scale, ~ 10 16 GeV [35]. In a variety of string theory models [35, 36], 
the PQ energy scale results in the range 10 16 GeV < f a < 10 18 GeV. From Fig. 1 we see that 
this range of f a values cannot be reconciled with axions as 100% CDM in Scenario I, while 
they can be in Scenario II provided Hi < 10 9 GeV. 

We studied the possibility for the axion to form the totality of the cold dark matter in the 
light of WMAP5 observations. There are two scenarios for this to happen. In Scenario I, in 
which the PQ symmetry is restored after inflation (f a < Hi/2ir), one needs Hi > 4.57 x 10 11 
GeV and f a = (7.27±0.25) x 10 10 GeV, corresponding to an axion mass m a = (85±3) fieV. In 
Scenario II, in which the PQ symmetry breaks during inflation and is not restored afterwards 
(f a > Hi/2n), one can have f a >Ax 10 8 GeV (m a < 15 meV) and Hi constrained as in 
Fig. 1 (see Eq. (49)). Moreover, in this region, the misalignment angle 0, must be chosen 
appropriately, as Eq. (50) and Fig. 2. Large values of the PQ symmetric breaking scale 
f a then correspond to small values of the initial misalignment angle d; t . Determining the 
PQ energy scale f a , or the axion mass m a , is thus tightly related to constraining inflation 
parameters and future extensions of the Standard Model. 
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